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Abstract. We investigate the analogy between circumstellar disks and the Taylor-Couette flow. Using the Reynolds 
similarity principle, the analogy results in a number of parameter-free predictions about stability of the disks, and 
their turbulent transport properties, provided the disk structure is available. We discuss how the latter can be 
deduced from interferometric observations of circumstellar material. We use the resulting disk structure to compute 
the molecular transport coefficients, including the effect of ionization by the central object. The resulting control 
parameter indicates that the disk is well into the turbulent regime. The analogy is also used to compute the effective 
accretion rate, as a function of the disk characteristic parameters (orbiting velocity, temperature and density). 
These values are in very good agreement with experimental, parameter-free predictions derived from the analogy. 
The turbulent viscosity is also computed and found to correspond to an a-parameter 2 x 10~* < a < 2 x 10~^. 
Predictions regarding fluctuations are also checked: luminosity fluctuations in disks do obey the same universal 
distribution as energy fluctuations observed in a laboratory turbulent flow. Radial velocity dispersion in the outer 
part of the disk is predicted to be of the order of 0.1 km/s, in agreement with available observations. All these 
issues provide a proof of the turbulent character of the circumstellar disks, as well as a parameter-free theoretical 
estimate of effective accretion rates. 
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1. Introduction 

Stars form by gravitational collapse of molecular clouds. 
During this process, proto-stars get surrounded and plau- 
sibly fed by the out-coming envelope/disk made of gas 
and dust, which can, under certain conditions, coagulate 
to form planetary embryos. There is little doubt that gas 
motions, usually considered as turbulent, play a major 
role. Turbulent motions enhance transport properties, 
thereby accelerate the evolution of the temperature 
and density in the envelope/disk. Also, turbulence may 
catalyze planet formation thanks to the trapping of dust 
particle inside large-scale vortices (Barge & Sommeria 
1995; Tanga et al. 1996; Chavanis 2000). As of now, the 
assertion that circumstellar disks are turbulent (what we 
shall refer to as the " turbulent hypothesis" ) has however 
never been properly checked. It mainly relies on the fact 
that the luminosity produced by the disk interacting with 
the central star is very large (see e.g. Hartmann et al. 
1998). In certain cases (FU Orionis-type systems), this 
luminosity is so high that it even supersedes the stellar 
component. The most widely accepted scenario so far to 
account for the abnormal luminosities of young stellar 



objects involves a magnetospheric accretion for classical 
T Tauri stars. In this case, the matter in the inner parts 
of the disk is coupled with the stellar magnetic field 
and falls onto the stars along the field lines at the free 
fall velocity (see e.g. GuUbring et al. 1998, MuzeroUe et 
al. 1998). This entails an accretion shock at the stellar 
surface in which almost all the visible and UV excess 
is produced (veiling the stellar lines). This scenario has 
shown many successes in interpreting spectral features 
of T Tauri stars like sodium and hydrogen line profiles 
(see e.g. MuzeroUe et al. 1998). In the case of Fu Ori 
stars, the currently advocated picture involves a wide 
boundary layer (see e.g. Popham et al. 1996). One of 
the main weaknesses of this scenario is its low predictive 
power since it all relies on an adjustable parameter (the 
accretion rate) which must be postulated a posteriori by 
comparison with observational data. 

Indeed, in this framework, the luminosity is directly 
related to the amount of energy released by the disk. 
The necessity for turbulence comes from the hypothesis 
that no laminar motions can produce the amount of en- 
ergy dissipated required to explain observed luminosities 
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(see e.g. Pringle 1981 and references therein). However, 
no attempt has ever been made to substantiate this claim 
in a quantitative manner. Questions we address here are: 
what are the luminosities produced by a laminar disk and 
by a turbulent disk ? and how do these compare with ob- 
servations ? These two questions are equally difficult to 
answer, but they hide different levels of difficulties. One 
is of theoretical nature: the physical processes at work in 
this disk/star interaction region are complex. A correct de- 
scription should include simultaneously the resolution of 
turbulence (with compressibility effects), radiative trans- 
fer (accounting for UV-irradiation by the star), magnetic 
processes, chemistry, the disk flaring, phase separation, 
time evolution, etc. The second difficulty is of observa- 
tional nature. At the present time, information has been 
gathered about the temperature, density and velocity dis- 
tributions in the outer parts of disks, at > 100 A.U. typi- 
cally, thanks to high resolution interferometry and clever 
data analysis (Guilloteau & Dutrey 1998). Unfortunately, 
basic parameters (mean free path, sound velocity and vis- 
cosity) connected with the gas dynamics and dissipation 
are still not known in the inner regions, and especially in 
the region where the disk and the star interact. 



Because of these difficulties, we choose to adopt radi- 
cally different approach than the classical model: instead 
of trying to build a fully "realistic" circumstellar disk, we 
use a simplified hydrodynamic model (" zero order model" ) 
and study in detail its physical properties. In the fu- 
ture, we will slowly increase its complexity (and reality!) 
by adding new ingredients like magnetic field, stratifica- 
tion, radiative transfer. Here, we show that our zero order 
model is analog to an incompressible rotating shear flow. 
It is therefore amenable to a simple laboratory prototype, 
the Taylor-Couette flow. From theoretical and experimen- 
tal studies of the properties of this prototype, one can then 
build general laws in circumstellar disks by a simple use of 
the Reynolds similarity principle. Taylor-Couette flow is 
a classical laboratory flow, and it has been the subject of 
many experiments. A recent review about stability proper- 
ties and transport properties for use in astrophysical flows 
has been made by DubruUc et al (2004b). As a result, they 
derive critical conditions for stability, and simple scaling 
laws for transport properties, including the influence of 
stratiflcation, magnetic field, boundary conditions and as- 
pect ratio. In the present paper, we apply these results to 
circumstellar disks and derive the expression of the char- 
acteristic parameters of the model as a function of astro- 
physical observables. We propose a procedure of quantita- 
tive estimate of the observable using observational results 
of Guilloteau & Dutrey (1998) and derive parameter-free 
predictions about turbulence and turbulent transport in 
circumstellar disks. These predictions are tested against 
observational data from T Tauri and FU Ori stars. 



2. Hydrodynamic model 

2.1. Basic ideas 

Observation of circumstellar disks suggests that they have 
sizes between 100 and 1000 astronomical units. In the se- 
quel, we will focus only on the part of the disk expected to 
behave like an incompressible fluid. An estimate of the im- 
portance of compressibility can be obtained via the Mach 
number, the ratio of the typical velocity to thermal veloc- 
ities. It is generally admitted that compressibility effects 
start playing a role when this number reaches values of 
unity. In the outer part of the disk, this ratio has been 
estimated by Guilloteau and Dutrey (1998) from CO line 
profiles. Its value is about 0.2-0.3. In the inner part of 
the disk (radius ranging from 1 to 30 astronomical units), 
we may use the disk structure inferred from the D/H ra- 
tio measured in the Solar System (Drouart et al, 1999, 
Hersant et al. 2001), which leads to a Mach number of the 
order of 0.05 to 0.1. These figures indicate that both the 
inner and the outer part can be treated as incompressible 
fluids. Closer to the star, the situation is less clear. On 
one hand, temperature tends to increase strongly, leading 
to an increase of the sound velocity and a decrease of the 
Mach number. On the other hand, as one gets closer to 
the boundary, one may expect larger typical velocities in- 
duced by larger velocity gradients, and thus increase of 
the Mach number. There are no direct observations sup- 
porting one scenario or the other. We shall then consider 
two scenarii: one, in which the Mach number Ma never 
exceeds unity. In this case, the whole disk is incompress- 
ible, and connects smoothly onto the star at the star ra- 
dius. The inner boundary is thus defined as r.i — r^. In a 
second scenario, the Mach number reaches unity at some 
"interaction radius" rjn, leading to an inner boundary at 
Ti = Hn- At this location, a shock appears, in which all 
velocities are suddenly decreased to very small values. In 
the shock, all the kinetic energy is transfered to the ther- 
mal energy, thereby producing a strong temperature in- 
crease (by a factor of 1 -f Ma^ ~ 2). This entails an in- 
crease of ionization and the matter gets more coupled to 
the stellar magnetic field. This second situation is consid- 
ered in magnetospheric accretion models (Hartmann et al, 
1998, Hartmann et al 2002), in which case is the Alfven 
radius; see Schatzman (1962,1989). Figure ^ summarizes 
these two possible configurations. From a hydrodynami- 
cal point of view, in the first situation the boundary is 
similar to free-slip boundary (with possible non-zero ve- 
locities in the direction tangential to the star boundary), 
while in the second situation, the interaction radius acts 
as a no-slip boundary (with all velocities becoming zero). 
This difference may refiect in the transport properties, 
see DubruUe et al. (2004b). In the sequel, the free-slip 
boundary condition will be referred to as smooth, while 
the no-slip boundary will be referred to as rough. In the 
laboratory experiments reviewed in DubruUe et al (2004b), 
the turbulent transport depends on the boundary condi- 
tions. Specifically, transport is enhanced (with respect to 
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any other boundary conditions) with boundary conditions 
of the rough or no-shp type. In the astrophysical case, it 
is not quite clear whether these two boundary conditions 
apply, or even whether different inner and outer boundary 
conditions result into an intermediate transport enhance- 
ment. We shall therefore devise observational tests using 
quantities independent of boundary conditions via a suit- 
able non-dimensionalisation. 



J L 

boundary 
/ layer 




b) indirect interaction 

Fig. 1. Two possible configurations considered in the 
present model: (a) the whole disk is incompressible and ex- 
tends onto the proto-star, and (b) the disk is incompress- 
ible until an " interaction radius" imposed for instance by 
a magnetic field. 



2.2. Basic equations 

In any case, the angular velocity il at the inner bound- 
ary is that of the star, namely ri(ri) — fJ*. For r > r^, 
the Mach number of the flow is less than one by construc- 
tion, i.e. pressure fluctuations vary over a time scale short 
compared with the dynamical time. In such a case, one 
can assume hydrostatic equilibrium in the vertical direc- 
tion, implying a decoupling of the vertical and horizontal 
structure. It is then convenient to describe the disk by its 
"horizontal equations" , obtained by averaging the original 
equations of motion in the vertical direction. The proce- 



dure is described e.g. in DubruUe (1992). It leads to: 

= M(5,S + 9^.) + (C + ^)a>ii) 



Here, is the horizontal gradient [dz — Q), u and p are 
the Favre average of the velocity and the pressure over the 
vertical direction, S is the surface density, /i and C are sur- 
face viscosity coefficients, G is the gravitational constant, 
M the mass of the star, and r the distance to the star in 
cylindrical coordinate, is a unit vector in the radial di- 
rection and H the vertical scale height. In the hydrostatic 
approximation, 



H ■ 



GM 



1/2 



(2) 



where Cg is the sound velocity. This expression is only 
valid when self-gravity can be neglected, namely when: 

Mdisk ^ H 

~M~ - R ^-^^ 
where Mdisk and M are the masses of the disk and the 
star, respectively (see e.g. Hure, 2000). 
In the opposite case, H will rather vary like the Jeans 
length in the vertical direction, as: 

These equations should be supplemented with an equa- 
tion for the surface energy E ~ Hc^, but we shall not need 
it in the sequel. 



2.3. Stationary axi-symmetric state 

The equation admits simple basic state, under the 
shape of stationary axi-symmetric solution. The mass con- 
servation then implies: 



or 



—dr{r'Eur) — 0, 
r 



(5) 



(6) 



where Mt is a constant, dimensionally equivalent to an 
accretion rate. Plugging this into the radial and azimuthal 
component of Q, we obtain two equations: 



Ml 



1 



rdr^ 



— ^ -drHp - 
r 



GM 



rdr^ 



S2 



The general solution of the second equation of Q is : 



(7) 



ug — A exp 



{P+l)dx/x] (8) 
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where A and B are constants and (3 — Mt/fi. Plugging 
this solution into the first equation of ^ then defines the 
general pressure. 

The basic state depends on three constants A, B and 
Mf, which must be specified through some sort of bound- 
ary conditions. In the case of astrophysical flows, the 
boundary conditions are not very well known and it is less 
easy to constrain the parameters. The condition that the 
rotation velocity of the disk matches the star velocity at 
the interaction radius only provides one relation between 
the three parameters: 

A^-^nVL^. (9) 

An additional constraint comes from the hypothesis that 
the disk is geometrically thin. This is consistent with our 
assumption that the Mach number is less than unity. In 
that case, as soon as there is no dramatic variation of the 
thermodynamic variables, radial pressure gradients and 
terms involving the radial velocity^ can be neglected in 
front of gravity, and the only way to satisfy the first equa- 
tion of (|7J) is to set Mt/fi = j3 = —3/2. In this case, the 
disk is almost Keplerian and obeys: 

3/i 



ue ^ ^^^ + ^{n,-nK{r,)), (10) 

where VLK{ri) is the angular Keplerian velocity at r = Vi. 
The corresponding solution is plotted in Fig. [3 It is made 
of a Keplerian disk in the outer part, with a continuous 
matching towards the star velocity at the interaction ra- 
dius. 

2.4. Comment about the disk outer radius 

In the expression we derived for the velocity, we did not 
need to specify anything about the condition at the outer 
boundary of the disk. One may then wonder what deter- 
mines this boundary, and whether it is relevant in spec- 
ifying the geometry of the problem. One way to answer 
this question is to note that stationary solutions of the 
shape H1U|I arc only possible provided the dynamical time 
scales are short with respect to the viscous time scale. In 
the vertical direction, the dynamical time scale to ensure 
hydrostatic equilibrium is H/cs ~ ^k^- In the horizontal 
plane, the two dynamical time scale are the radial time 
scale r/ur ~ r'^'S/fi and the orbital time scale fl'^ ^ ^k^- 
The radial time scale is comparable with the viscous time 
scale. So the condition for stationarity is that the orbital 
time scale is less than the viscous time scale, resulting in 
r < To, with To solution of the equation: 

^ The viscosity and the advection terms are then negligible 
in front of the radial pressure gradient by a factor of the order 
of Ma— and Ma^, respectively. 




r (A.U.) 



Fig. 2. Velocity profile in a circumstellar disk in the vis- 
cous regime, with keplerian velocity in the outer part, and 
smooth matching onto the central object at the interac- 
tion radius. For this example, the mass of the central star 
has been taken as a solar mass. 



This radius defines the outer geometrical limit within 
which stationary solutions can reasonably exist. If in this 
formula, one considers the ordinary viscosity, then one 
typically obtains Tq much larger than the observed disk 
outer radii. On the other hand, one may argue that as 
soon as the disk is turbulent, the molecular viscosity be- 
comes irrelevant, one must consider a kind of "turbulent 
viscosity" in this formula. In this case, using the formula 
we derive in Section 3.3, we find that is of the order of 
the disk scale height. In geometrically thin disks, it is not 
quite clear whether this limit really exists or not. Neither 
limits really match the observed disk radii. However, the 
sharp observed edge of disks remains inconsistent with 
stationary models. Stationary solutions, due to constant 
accretion rate in radius, are indeed in essence radially in- 
finite. This suggests that the disk outer radii may still be 
linked with some unstationary effects. This is the subject 
of a forthcoming paper (Mayer et al. 2004, in preparation). 

2.5. Ttie incompressible analog 

Astrophysical disks are (weakly) compressible and radi- 
ally stratified. It is however possible to build an incom- 
pressible analog of them, using clever boundary condi- 
tions. This remark is at the heart of the laboratory pro- 
totype. Consider indeed an incompressible, unstratified 
fiuid, enclosed within a domain with cylindrical symmetry, 
bounded by inner and outer radii and Tq. Its equation 
of motions are given by the Navier-Stokes equations: 

dtu + u-Vu = — Vp + i/Au, 
P 

V • u = 0. (12) 



Hersant et al.: Turbulence in circumstellar disks 



5 



where p is the density, u is the velocity, v the molecular 
viscosity and p is the pressure. If we assume hydrostatic 
equilibrium in the vertical direction, we get: 

9^P^0, (13) 

so that p is a function of r only. In that case, equation (|12|l 
admits simple basic state under the shape of stationary 
solutions, with axial and translation symmetry along the 
disk rotation axis (the velocity only depends on r). The 
incompressibility condition then implies: 



or 



— dr{rUr) ~ 0, 

r 

K 



(14) 



(15) 



where K is a, constant, to be constrained later. Plugging 
this into the radial and azimuthal component of H12|l . we 
obtain two equations: 



r3 



r p 

K, 

2' 



r 



:dr{rue) = {^^^r^^ 



(16) 



The second equation of (|16() is homogeneous in r. It only 
admits two power law solutions, with exponent —1 and 
1 + K/v, so that the general solution is: 

(17) 



ue = Ar^+l^ + - 
r 

where A and B are constants and (i = K/v. Plugging this 
solution into the first equation of l|16|l then defines the 
pressure. 

The basic state depends on three constants A, B and 
which must be specified through boundary conditions. 
In laboratory fiows, these conditions are usually well de- 
fined and allow for a simple determination of the constants 
once the rotation velocities at the inner and outer bound- 
aries are known (Bahl, 1970): 



A 



B = 



I _ ^/3+2 
^2 



1 



13+2 



(18) 



where 77 = ri/vo is the radius ratio, flo and fli are the 
angular velocity at outer and inner radii and (3 = Rj. ~ 
K/v = Ur{ri)ri/v is the radial Reynolds number, based on 
the radial velocity through the wall of the inner cylinder. 
Note that it is positive for motions outwards from the axis 
of rotation. 

Comparing H10() with H18|l and (|17() . it is possible to 
see that the "laboratory" analog of Keplerian fiow is such 
that: 

3 



n,. 



2' 



IGM 



n/ro- 



(19) 



This shows that this basic state describes as well labora- 
tory incompressible flows, with rigid boundaries and with- 
out gravity, in which angular momentum distribution may 
be imposed by boundary conditions, and astrophysical 
flows, without rigid boundaries, in which angular momen- 
tum distribution is imposed by gravity. In other words, 
building a prototype of astrophysical disks (within the 
approximations described above) requires a (laboratory) 
flow with equivalent angular momentum distribution, and 
equivalent control parameters. We now derive these con- 
trol parameters, in order to apply the Reynolds similarity 
principle. 

2.6. Control parameters 

The shape of the basic state allows for the determination 
of the control parameters of the flow. These parameters are 
essential in the comparison with the laboratory prototype 
since the Reynolds similarity principle states that the as- 
trophysical disk will behave like the laboratory prototype 
with same control parameters. These control parameters 
are 

the global Reynolds number: 

S{ro ~ n? 



Re 

the rotation number: 
the curvature number 



Rn — 



2Vt 



Rc = 

ro - Ti 

the local radial Reynolds number: 



Rr 



the aspect ratio: 



F = 



H 



(20) 



(21) 



(22) 



(23) 



(24) 



Here, fl, S and r are characteristic angular velocity, shear 
and radius. Adopting the convention of Dubrulle et al 
(2004b), we find: 

\ 1/2 
To I r„ 



— 1 / T \ ^ /'^ — 



(25) 



while f is fixed through the condition Q — ug(f)/f. A 
simplification occurs in two limiting cases, relevant to as- 
trophysical disk: Ti ^ To or ri 17* flxiri). In 
both cases, we have: 



2/3 1/3 

s - -^nKif). 



(26) 
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The control parameters then siniphfy into: 



Re = 



Rn = 



Rc = 



Rr — — X 



4 
3' 

2/3 1/3 

r.- rn 



where T is the disk temperature, and n the number den- 
sity. 

When the gas is weakly ionized (xg <C 1), the transport 
coefficient must be multiplied by a factor (Lang, 1980): 



4 X 10 



-12 



y2 



(31) 



3 
2' 
H 



This correction is valid as long as Vionjvneu < 1- This 
sets a limit on ionization fraction, below which the gas 
viscosity takes the neutral value: 



2/3 1/3' 



(27) 



jif = 4 X lO^^^T^ 



(32) 



If the disk is stratified or magnetized, other control 
parameters appear, like the Prandtl number Pr = v/k 
and the magnetic Prandtl number Pm — ^ov/rj where k 
and rj are the heat diffusivity and the magnetic resistivity, 
and /ig is the vacuum permeability. 

2. 7. Molecular transport processes 

The computation of the control parameters requires an 
estimate of molecular transport coefficients. These coeffi- 
cients depend on the ionization state of the gas. There are 
two sources of ionization. Thermal ionization is efficient in 
the inner part of the disk. The corresponding ionization 
fraction can be written as (Fromang et al, 2003): 



The Prandtl number in this case is equal to 10 (Lang, 
1980). 

The resistivity of an ionized gas can be written as the 
sum of the resistivity induced by electron- neutral collisions 
and electron-ion collisions: 



V = Ven+ Ve 



where (Lang, 1980) 



= 10" 



'■-T^/'^ohm - 



nln A 
4 X lO'^^rj- ohm - 



y3/2 



(33) 

(34) 
(35) 



„th 



6x10"^ 



T 



\ 3/4 

1000/f ) 



2.4 X lO^^cm" 



1/2 



exp(- 



-251^^1^/:^): 



= 1.3 X 10 



4 7^3/2 



is the Coulomb logarithm. 



(28) 

where T and n are the temperature and the number den- 
sity of neutral species (hydrogen mainly), respectively. For 
temperature lower than 10'^ K (typically r > 1 A.U.), ther- 
mal ionization is negligible. However, X-Ray illumination 
from the central star (or the magnetospheric accretion 
flow) may induce a weak ionization in some part of the 
disk (typically away from the mid-plane) (Feigelsson and 
Montmerle, 1999). A recent theoretical study has recently 
been performed by Igea and Glassgold (1999). They found 
that at a given radius from the source, the ionization frac- 
tion is a universal function of the vertical column density 
N± , independent of the structural details of the disk. The 
role of cosmic rays in the disk ionization is still a matter of 
debate (Sano et al. 2000) and will not be discussed here. 
For a typical young stellar object, Igea and Glassgold's 
result can be approximated by: 



2.8. Physical parameters 



n 



_rpl/4^-l/2^-0.002(rAu-l) 



37.1/4 -1/2 -0.002(rA!7-l) 



where Ue is the ion number density, and rAu is the distance 
from the central star, in astronomical units. 

When the gas is neutral, the viscosity and heat diffu- 
sivity are given by (Lang, 1980): 



= K = 3 X 10 



7.1/2 



Various physical parameters are required to estimate the 
control parameter. 

Parameters associated with the disk are Tq, r^, F and 
v. The disk inner radius depends on whether the disk/star 
interaction is direct or indirect. In the first case, = r*. In 
the second case, r.^ may no exceed the corotation radius, 
at which the disk velocity matches the star velocity. In 
the sequel, we shall consider variations of in between 
these two limits. The disk outer radius To must be specified 
through the implicit relation (|ll|l . Direct observation for 
disk suggest that the disk size is of the order of r^i = 1000 
A.U. for disk around T Tauri and even smaller ro — 50 
A.U. for disk around FU ORI (Kenyon, 1999). Clearly, ro 
is thus the maximum size Tq can achieve. For practical 
reasons, we defer its discussion in the next section, after 
computation of the temperature and density profile. 
N± > 10 ciparameters relative to the proto-star have been mea- 
sured in some T Tauri and FU Ori stars. Table ^ gives 
< 10 cm aisai(:§Sljb of stars we shall use in the following. It is par- 
ticularly convenient and illustrative to scale all quantities 
in the problem with respect to values defined at the dis- 
tance of f = 0.33 A.U., which is the characteristic radius 
corresponding to a disk with = 10"'^^ cm and Tq — 1000 
A.U., the two extreme limits for and Tq. We also use, as 
a reference, a rotation period of the star of 8 days (typical 



(30) T Tauri star), leading to 1^* = 9 x lO'^ s" 
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Temperature and number density in circumstellar 
disks are not known due to the lack of spatial resolu- 
tion. However, their magnitude can be deduced by short 
radii extrapolation of measurements made on the outer 
disk. The inversion method of Guilloteau & Dutrey (1998), 
based on ^-square fitting of CO interferometric maps, 
yields the temperature and the density profiles at r > 
50 — 100 A.U.. For instance, for the disk around DM Tau 
{M ~ O.5M0), their method predicts 



r~ 30 



io« 



(100 Au) 

r 



-2.75 



1 



100 AU. 

r \ -1-5 



(100 Au) 



K 



g cm 



.100 AU. 

where S is the surface density of the disk (including gas 
and dust). Error bars on the measurements are rather 
large, and could amount to a possible variation by a fac- 
tor 5 to 10. At r = 1 A.U., one finds n ~ 3 x 10^^ 
cm~^ and E ~ 1000 gcm~^, in agreement with values ob- 
tained by modeling the deuterium enrichment in the Solar 
System (Drouart et al. 1999; Hersant, Gautier & Hure 
2001). At the reference radius of 0.33 A.U., the density is 
n ~ 7x 10^** cm^'^, the surface density is S ~ 5275 gcm^^, 
the height is H — 4.6 x 10^^ cm and the temperature is 
~ 925 K. In a more recent analysis of the disk around BP 
Tau, Dutrey at al. (2003) found densities and temperature 
corresponding to a value of n ~ 3 x lO^"' cm^'^, S ~ 992 
g cm-2, H = 10" cm and T = 289 K. The difference 
between these figures and the figures of DM Tau provide 
an illustration of the error bars associated with our "typ- 
ical values", since the disk around BP Tau seems much 
smaller, and correspond to a more evolved stage than the 
disk around DM Tau. 

2.8.1. Ionization state 

It is interesting to study the ionization state of the disk 
with temperature and density observed in DM Tau. The 
ionization fraction is plotted as a function of radius in Fig. 
^ for the thermal and X-ray contribution. One sees that 
the thermal contribution dominates at r < 1 A.U., while 
the X-Ray contribution becomes important at larger radii. 
However, comparing with the limiting ionization state eq. 
H32(l . one sees that only the outer part of the disk r > 
100 A.U. is sufficiently ionized to infiuence the molecular 
viscosity. 

2.9. Regimes 

The previous scaling allows for an estimate of the disk con- 
trol parameter, and, therefore, for an identification of the 
possible regimes. A difficulty with respect to the labora- 
tory experiment is that in disks, the transport coefficients 
vary over across the disk due to the radial stratification. 
To define the control parameter, one must pick up a typ- 
ical value. In this context, it is logical to consider their 
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Fig. 3. Midplane ionization fraction due to thermal con- 
tribution (dot-dashed line) and X-ray contribution (dot- 
dotted line). The plain line is the limiting ionization frac- 
tion, below which the ionization does not influence molec- 
ular transport. The shaded area is the region where ion- 
ization has to be taken into account in the computation 
of the viscosity. 

value at r = f , since both the typical shear and radii at 
this location have been used. Using the values given in the 
previous section, we then obtain: 



i?e = 2 X 10 



25 



M 



1/2 



T 



10" cm 



7 X 10" cm-3 J V 930 K 

-4/3 / X -2/3 



-1/2 



103 A.U. 



Pr 
Pm 

r 



4 

~3' 
3 

"2' 
1, 

2 X 10" 
0.94. 



(36) 



The value of the rotation number indicates that the 
flow is anti-cyclonic and belongs to the "globally sub- 
critical" class defined in DubruUc et al (2004b) . The radial 
Reynolds number is negative, indicating an inward radial 
circulation. Its value is close to unity. So its influence on 
transport properties can be neglected as a first approxi- 
mation, see DubruUe et al (2004b). The curvature number 
depends on the interpretation of the viscous time scale (see 
Section 2.9.2). In the case where the interpretation is done 
with the molecular viscosity, one finds Rq ~ 0.0004, that 
is disks are in the wide gap limit. In the other limiting 
case where the viscous time scale is computed using the 
turbulent viscosity, one finds Rc = 1 — F, that is disks are 
in the small gap limit. The relevant Reynolds parameter 
to be considered in studying transport properties will be 
ReRl; = Sf'^/v, see DubruUe et al (2004b). Finally, the 
aspect ratio is less than one. From the review of DubruUe 
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et al, we infer that an additional correction must be 
included in the definition of the relevant Reynolds param- 
eter, which becomes: 



Re. 



phys 



ReiRcT) 



-1/2 



n 



7 X 1014 

To 

103 A.U. 



T-3 



-1/2 



T 



930 K 



This is the expression one would naturally derive by con- 
sidering the "smallest" length scale in the problem, see 
e.g. Longaretti (2003). 



3. Predictions about the structure of circumstellar 
disks 

3.1. Stability: the laminar/turbulent transition 

The stability properties of circumstellar disks can be found 
by comparing the physical Reynolds number Repays with 
critical Reynolds numbers derived in laboratory experi- 
ments, in the anti-cyclonic non-linear regime. These mea- 
surements are summarized in DubruUe et al (2004b). 
Disregarding any body forces, one finds a critical Reynolds 
number of the order of 2300, well below the disk value. 
Taking into account the possible stable vertical stratifi- 
cation observed e.g. in DM Tau (Dartois et al, 2003) ^, 
one obtains a slightly larger value of the order of 4000 
(DubruUe et al, 2004a). The presence of a vertical mag- 
netic field may increase the critical Reynolds number, due 
to the low magnetic Prandtl number prevailing in disks. 
Using the scaling of Willis and Barenghi (2002), one finds 
a critical Reynolds number Rcc ~ lOQ/Pm ~ lO^". This 
is still well below the observed Reynolds number. These 
number concur to conclude that the disk is turbulent. 

However, due to the huge variation of the transport 
coefficients across the disk, one may wonder how strong 
this conclusion is. A way to answer this question is to 
see how locally in the disk, the stability criterion are sat- 
isfied using "local" non-dimensionalized parameter, built 
by replacing f by r the distance to the central object. The 
result of this procedure is plotted on Fig 01 One sees that 
at any radius, such effective local Reynolds number is well 
above any critical Reynolds number due to body forces. 
This strengthens our conclusion. 

3.2. Mean energy dissipation and accretion rate 
3.2.1. Definition 

The quantitative comparison between experimental mea- 
surements and energy dissipated in circumstellar disks 
first requires a relation linking the torque and the disk 



^ This may be due to disk illumination by the central star 
(D'Alessio et al, 1998) 
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Fig. 4. Physical local Reynolds number in circumstellar 
disks as a function of radius (dotted- line with symbols). 
The dashed-dotted line and the full line are the critical 
Reynolds number deduced from laboratory experiments, 
see DubruUe et al (2004b). 

luminosity C. The total power dissipated in a Taylor- 
Couette experiment with same control parameter as in 
a keplerian disk is 



SG 



G 



(38) 



where S is the disk surface density, G the non-dimensional 
torque and D the viscosity. In a stationary disk, this power 
is dissipated under the form of heat, and thus coincides 
with the disk luminosity, that is 



(39) 



In practice, observed luminosities are often expressed as 
a function of an "effective mass accretion rate", namely 
(Hartmann et al. 1998) 



M 



GM 



(40) 



From a theoretical point of view, the detailed computation 
of the accretion luminosity is not straightforward since it 
depends on the boundary condition at the interaction ra- 
dius. For comparison with experimental data, we therefore 
consider the quantity: 



G 



Re^ 



M 

Mo' 



where Mq is an effective accretion rate given by 



H 



Mo = Sr,ff7 — 



(41) 



(42) 



The quantity G/Rel^y^ (the non-dimensional energy dis- 
sipation) , includes all the boundary condition dependence 
and only depends on the Reynolds number. 
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3.2.2. Prediction using laboratory experiments 



The result of DubruUe et al (2004b) lead to an analytical 
prediction for the function M/Mq = f{Rn)Gi/Rep^y^, 
where f{Rn) is a function parameterizing the influence of 
rotation, and Gi is the torque in situation when only the 
inner cylinder is rotating. From the experiments, we infer 
^ 0.1 if the flow is turbulent, and 1 is the flow 
is laminar. Taking the wide gap limit <C Tq, R^Re = 
Rsphys in the formulae of Dubrulle et al (2004b), we obtain 
three possible regimes: 



— in the laminar regime, for Rephys < 2300 



M 

Wq 



27r 



R^phys 

for smooth boundary conditions and Rephys > 2300 



(43) 



M 



oth 



0.06 [ln(3 



10 '^Re^phys) 



I -3/2 



(44) 



for rough boundary conditions and Repays > 2300 



M 

Mr 



I rough 



= 1.9 X 10" 



(45) 



In the other limit — s- Tq, the turbulent value take 
similar expression, but must be multiplied by a factor 
1/3(1 — r)'^/^. In astrophysical disks, the boundary con- 
ditions are not known a priori. Moreover, given the huge 
physical difference between the inner part and the outer 
part, it is unlikely that the boundary condition at the in- 
ner and outer part coincide, so that we are probably more 
in a state of " mixed" boundary conditions studied exper- 
imentally by Van den Berg et al (2003). In that case, the 
energy dissipation is found to vary in between the two lim- 
its set by respectively the "pure" smooth type (|44|l and 
the pure "rough" type (03, see Dubrulle et al (2004b). 
We shall therefore adopt these formulae as a lower and 
upper limit of the energy dissipation in disks. 



3.2.3. Test against observational data 

We find the following scaling for Mq: 



Mo ~ 3 X 10" 



5300 gem 



M \ ''"/ \ 
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Fig. 5. Comparison between the non-dimensional energy 
dissipation predicted from laboratory measurements in the 
wide gap limit (dotted lines) or in the small gap limit 
(plain lines) and observed in circumstellar disks (sym- 
bols) , as a function of the Reynolds number Rcphys ■ For an 
easier comparison, the mean energy dissipation has been 
translated into the non-dimensional accretion rate M /Mq 
(computed using Ea. (|m and observationally determined 
parameters reported in Tab.QJ. The symbols ffl report the 
value using ri — rcoro, which provides an upper bound of 
the energy dissipation. The circles report the value using 
ri — r, , which provides the lower bound of the energy 
dissipation. All the quantities have been computed using 
temperature and density estimated for the DM Tau sys- 
tem using the results of Guilloteau and Dutrey (1998), so 
there is no adjustable parameter in this plot. 



We are aware that these values are probably uncertain by 
a factor of 10 or even more. In the next decade, the results 
expected with ALMA will largely shorten the error bars. 

At i?e = 3 X 10^3, our model (Eqs. and (jlHl 

predicts that 0.0002 < M/Ma < 0.019, resulting in 
6 X 10"*^ < M < 6 X IO-HIqjt-'^ for T Tauri stars. This 
is in good agreement with the observed values ranging 
from M = 10-1° to lO-^Mgyr-i (Hartmann et al. 1998). 
Disks around FU Ori are characterized by a smaller disk 



radius, leading to higher values of Al according to (|46|) . 
by a factor 5. This is not quite enough to reach values of 
up to lO-^M© yr-i, associated to disks around FU Ori 
stars (Kenyon 1994). Such values could be obtained if the 
typical disk density is higher in disk around FU Ori, re- 
sulting in more massive disks. This is plausible, since FU 
Ori are younger than T- Tauri stars. 

A graphical representation of this discussion can be 
obtained by plotting the computed M / Mq as a function 
of Rephys using the values listed in Table ^ as input pa- 
rameters. To remove the problem with our ignorance of the 
MQC4&Kctual value of r^, we have used the relation r* < < rcoro 
and computed the corresponding M / Mq and Rephys- The 
actual dissipation somehow lies in between the two corre- 
sponding estimates. These estimates are plotted in Figure 
|S1 For comparison, we have added the theoretical predic- 
tions (eq. (|44|l and l|45|) giving the minimum and the max- 
imum expected values in the turbulent case, as well as the 
laminar value. This last value is very much lower than the 
turbulent values, and is never even nearly approached by 
any stars we considered. This may be seen as a proof of 
the turbulent character of all these disks. 



10 



Hersant et al.: Turbulence in circumstellar disks 





input parameters 


output parameters 




mass 


radius 


coro. radius 


disk outer edge 


accr. rate 


eff. rates 


Reynolds num. 


star 


Ml Mr,. 


' * / J- "-t^ 


' COl O / ''J 


ro/1000 A.U. 


log M 


lop- Mn(r.A 


lop- Mn(r +^ 

lug iv±{j yi corot J 


log 


lug 1 tccoro 


AA Tau 


0.53 


1.74 


13.81 




-8.48 


-4.39 


-3.67 


13.53 


12.55 


BP Tau 


1.32 


1.99 


17.79 




-7.54 


-2.87 


-2.10 


13.27 


12.24 


CY Tau 


0.42 


1.63 


10.05 




-8.12 


-4.17 


-3.54 


13.61 


12.74 


DE Tau 


0.25 


2.45 


10.05 




-7.58 


-4.00 


-3.52 


13.54 


12.85 


DF Tau 


0.27 


3.37 


11.29 




-6.91 


-3.31 


-2.90 


13.39 


12.78 


DK Tau 


0.43 


2.49 


10.05 




-7.42 




O.Ul 


13.42 


1 9 7'^ 


DN Tau 


0.38 


2.09 


10.03 




-8.46 


-4.60 


-4.06 


13.52 


12.76 


DO Tau 


0.37 


2.25 


10.05 




-6.84 


-3.00 


-2.48 


13.50 


12.77 


DQ Tau 


0.44 


1.79 


10.05 




-9.40 




-4 


1 ^ ^fi 

±o.ou 


19 7^ 


DS Tau 


0.87 


1.36 


10.05 




-7.89 


-3.45 


-2.76 


13.53 


12.58 


GG Tau 


0.44 


2.31 


10.05 




-7.76 


O.O-L 


o.ou 


1 4'i 


19 7^ 


GI Tau 


0.71 


1.48 


10.05 




-8.02 


79 


— o.uu 


±o.oo 


1 9 (S^ 


GK Tau 


0.46 


2.15 


10.05 




-8.19 


-4.21 


-3.67 


13.47 


12.72 


GM Aur 


0.52 


1.78 


10.05 




-8.02 


-3.95 


-3.34 


13.52 


12.69 


HN Tau 


0.81 


0.76 


10.05 




-8.89 


-4.45 


~o.oo 


lO.OU 


1 9 (sn 


IP Tau 


0.52 


1.44 


10.05 




-9.10 


-5.00 


-4.33 


13.62 


12.69 


UY Tau 


0.42 


2.60 


10.05 




-7.18 


97 


-9 sn 


13.41 


12.74 


CI Tau 


0.5 


1.87 


10.05 




-7.19 


-3.14 


-2.56 


13.51 


12.70 


CX Tau 


0.33 


1.63 


10.05 




-8.97 


-5.18 


-4.55 


13.66 


12.79 


CZ Tau 


0.41 


1.19 


10.05 




-9.35 




-4 fi^ 


1 7^^ 

lO. 1 o 


12.74 


DM Tau 


0.43 


1.39 


10.05 




-7.95 


-"i Q7 


9Q 


1 fi7 


19 7^ 


DD Tau 


0.42 


1.44 


10.05 




-8.39 


-4.43 


-3.76 


13.66 


12.74 


DH Tau 


0.38 


1.67 


11.33 




-8.30 


-4.42 


-3.76 


13.62 


12.71 


DI Tau 


0.43 


1.71 


12.56 




-8.75 


4 7Q 


-4.01 


1 '\R 

lO.OO 


12.64 


DP Tau 


0.46 


1.44 


10.05 




-7.88 


-3.86 


-3.19 


13.64 


12.72 


FM Tau 


0.58 


1.17 


10.05 




-8.45 


-4.26 


-3.52 


13.68 


12.67 


FO Tau 


0.33 


1.59 


10.05 




-7.50 


-3.71 


-3.07 


13.67 


12.79 


FQ Tau 


0.35 


1.42 


10.05 




-6.45 


-2.61 


-1.93 


13.71 


12.78 


FS Tau 


0.46 


1.25 


10.05 




-8.09 


-4.06 


-3.34 


13.70 


12.72 


FV Tau 


0.71 


1.87 


10.05 




-6.23 


1 Q'^ 


-1 ^7 


13.44 




FX Tau 


0.34 


1.94 


10.05 




-8.65 


-4.86 


-4.28 


13.58 


12.79 


FY Tau 


0.50 


1.87 


10.05 




-7.41 


o.ou 


-9 7S 


1 ^ ^1 


1 9 7n 


GH Tau 


0.29 


1.90 


10.05 




-7.92 


-4.23 


-3.65 


13.62 


12.82 


GO Tau 


0.50 


1.40 


10.05 




-7.93 


-3.86 


-3.17 


13.64 


12.70 


Haro 6-37 


0.60 


1.90 


10.05 




-7.00 


-2.83 


-2.26 


13.46 


12.66 


HO Tau 


0.56 


0.94 


10.05 




-8.86 


-4.68 


-3.86 


13.79 


12.68 


IQ Tau 


0.35 


2.01 


10.05 




-7.55 


-3.74 


-3.18 


13.56 


12.78 


LkCa 15 


0.81 


1.53 


10.05 




-8.87 


-4.49 


-3.84 


13.49 


12.60 


Lk Ha 332/Gl 


0.29 


2.36 


10.05 




-6.60 


-2.93 


-2.43 


13.53 


12.82 


V955 Tau 


0.44 


2.34 


10.05 




-7.02 


-3.07 


-2.57 


13.44 


12.73 


V1057 Cygni 


0.50 


1.60 


13.53 


0.005 


-4.00 


-0.86 


-0.12 


14.73 


13.72 


FU Ori 


0.70 


1.20 


15.83 


0.005 


-3.70 


-0.31 


-0.58 


14.78 


13.58 



Table 1. Observational parameters for T Tauri stars (from Bouvier 1990, Hartmann et al. 1998) and for FU Ori stars 

(from Popham ct al. 1996) considered in this study (left) and disk physical parameter (right). The computation of the 
corotation radius requires the knowledge of the star rotation velocity. In case this last quantity is not available, the 
corotation radius has been set to 10.05, the solar value. Lower and upper bound on Mq and Re have been computed 
using either the star radius or the corotation radius for r^. Accretion rates are in M0yr~^. 



We also see that energy dissipation for disk around T 
Tauri has a tendency to cluster in between the minimal 
and maximal values allowed by the theoretical predictions. 
The relative position of the cluster of point is slightly bet- 
ter in the case where is computed with r corot, which 
may be an indication that is actually closer to the coro- 
tation radius than to the star radius. However, given the 
error bars stressed above, this is probably not enough to 



conclude that disk around T Tauri stars are connected 
through a magneto-sphere, rather than through a bound- 
ary layer. In the case of FU Ori, the points are clearly 
above the maximum allowed by our choice of parameters. 
The discrepancy is slightly lower for the case when = r* , 
a choice which should probably be favored by the possi- 
ble signatures of boundary layer in these objects (Kenyon, 
1994). In that case, an increase of the surface density by 



Hersant et al.: Turbulence in circumstellar disks 



11 



a factor 10 with respect to our values will be enough to 
solve the discrepancy with theory. 

Our estimate neglects the influence of the magnetic 
field. Laboratory experiments using liquid metals have 
proved that this can potentially change the intensity of 
the transport with respect to the pure hydro-dynamical 
case. However, no experiment has been performed so far, 
to study the magnetic influence in regimes relevant to as- 
trophysical disks. 

3.3. Turbulent viscosity 

The turbulent viscosity in circumstellar disks can be pre- 
dicted by comparison with laboratory measurements, see 
DubruUe et al (2004b). We find: 



1 Tia 



G 



2tt 



Re 



-\S\H^ 



(47) 



phys 



where r — Y,S is proportional to the angular momen- 
tum, and the index lam means laminar value. The ratio 
Tiam/T is & function describing the radial variation of the 
turbulent transport. In an incompressible laboratory fiow 
with constant density, this function is just the ratio of 
the laminar shear profile to the turbulent shear profile: 
the turbulence regulates itself through a modification of 
the velocity profile. In keplerian disks, the shear profile 
is fixed (through the gravitational force), and the reg- 
ulation operates through the density. This function has 
been measured in a number of laboratory experiments. 
At large Reynolds number, it seems to approach a con- 
stant value of 4 predicted by Busse (1970) using argument 
of maximal momentum transport. With \S\ — l.b^l and 
G/Re^i^yg — M / Mo, this defines a resulting typical tur- 
bulent viscosity as: 



draw interesting information from the luminosity fluctua- 
tions which refiect the dynamics of the underlying turbu- 
lent fiow. In laboratory experiments with smooth bound- 
ary conditions, turbulent fluctuations are observed to fol- 
low an universal (i.e. Reynolds number-independent), log- 
normal distribution (Lathrop, Fineberg & Swinney 1992) 
with variance A = 0.042. The universal distribution oc- 
curs for variable normalized by their mean. Energy dis- 
sipation is proportional to the wall shear stress squared. 
Since the functional shape of the log-normal distribution is 
unchanged by squaring, distribution of energy dissipation 
should also be log-normal. To check this prediction, we 
have computed the distribution of the luminosity fluctua- 
tions observed from the disk around BP Tau and from the 
disk around V1057 Cygni. The results are shown in Fig. 
El and [3 One sees that the fluctuations in the disk around 
V1057 Cyg are very well fitted by a log-normal distribu- 
tion, with a variance similar to that of laboratory experi- 
ments. In the case of BP Tau, however, the comparison is 
not as good. This difference between the two systems may 
be traced to different boundary conditions. If we accept 
that disk around T Tauri stars follow the magnctospheric 
accretion scenario, while the disk around FU Ori is con- 
nected to the star through a boundary layer, it may not be 
surprising that only the disk around FU Ori follow the lab- 
oratory, smooth boundary condition distribution. Since we 
do not have any measurements for rough boundary condi- 
tions, we cannot say whether the discrepancy comes from 
the different boundary conditions, or from the presence 
of other physical effects, like accretion shock, or magnetic 
field. 



Mo 



(48) 



This turbulent viscosity takes the shape of an a viscos- 
ity, proposed by Shakura and Sunyaev (1973). The corre- 
sponding a coefficient is here a function of the Reynolds 
number of the circumstellar disk (through eqs l|44(l and 
At Rcphys = 3x 10^^, its value is typically 2 x 10""^ < 
a < 2 X 10~^. This range is in good agreement with the 
range of values inferred from the D /H ratio in the Solar 
System (Drouart et al, 1999; Hersant et al, 2001). More 
generally, this range is compatible with disk lifetime and 
values usually adopted in theories. Using this expression 
for the turbulent viscosity, the viscous timescale writes : 




3 M 



(49) 



3.4. Energy fluctuations 



Up to now, we have considered only the mean energy dis- 
sipation and its luminous counterpart, but we can also 



Fig. 6. Distribution of luminosity ffuctuations observed 
disk around BP Tau (symbols) compared with a log- 
normal distribution of various variance A (plain line) . The 
value of A is indicated aside each line. 
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Fig. 7. Distribution of luminosity fluctuations observed 
disk around V1057 Cyg (symbols) compared with a log- 
normal distribution of variance A = 0.03 (plain line). 

3.5. Velocity fluctuations 

Laboratory measurements provide interesting clues about 
the intensity of velocity fluctuations. Since such fluctua- 
tions may be potentially observable in disks using non- 
thermal line widening, they may be used as additional 
constraints or observational test of the analogy between 
laboratory flows and circumstellar disks. From results of 
DubruUe et al (2004b), it appears that azimuthal veloc- 
ity fluctuations should be proportional to the mean az- 
imuthal velocity, with a proportionality factor depending 
weakly on the Reynolds number, like O.Q3(Re/ Rcc)"'^'^^^ . 
With Rcphys = 3 X 10^^ and Rcc = 10^ (see Section 3.1), 
the factor is of the order of 0.01. Using Eq. (|10(l . we can 
compute the azimuthal velocity dispersion for a typical 
circumstellar disk around a T Tauri. The azimuthal ve- 
locity dispersion decreases from about 0.6 km/s in the 
inner part, to 0.03 km/s in the outer part, at 100 A.U. 
The total velocity dispersion depends on the anisotropy of 
the turbulence. In the laboratory experiment, the radial 
relative velocity dispersion is observed to be about twice 
the azimuthal velocity dispersion. There was no measure 
of the vertical velocity dispersion, but it can be expected 
to be much smaller than the horizontal dispersion due to 
the rotation-induced anisotropy (DubruUe & Valdettaro 
1992). 

Velocity dispersion in disks have been measured by 
Guilloteau & Dutrey (1998) at r > 100 A.U. They obtain a 
value of the order of 0.1 km/s, which would correspond to 
a value of about 0.05 km/s for the azimuthal component. 
This is close to the values found from comparison with 
laboratory flows. 

4. Summary 

In this paper, we have derived and studied the analogy 
between circumstellar disks and the Taylor-Couette flow. 



This analogy results in a number of parameter-free pre- 
dictions about stability of the disks, and their turbulent 
transport properties, provided an estimate of the disk 
structure is available. We have proposed to get this es- 
timate from interferometric observations of circumstellar 
disks, and checked that the energy dissipation, the turbu- 
lent transport, and the fluctuations in circumstellar disks 
all follow behavior compatible with the prediction from 
the analogy. This check can first be used as a clear proof 
of the turbulent character of circumstellar disks. A second 
interesting application would be to build from this anal- 
ogy a parameter free model of circumstellar disks. In this 
respect, the proportionality between the turbulent viscos- 
ity and the so called "accretion rate" (a quantity easily 
accessible to observation) is very interesting because it 
opens the possibility to infer the disk structure from the 
observation of its luminosity. For this, a model has to be 
built linking the turbulent transport and the disk struc- 
ture. This is the subject of ongoing work. 

We note finally that our model could also possibly ap- 
ply to other type of disks (e.g. around black holes, or in 
close binaries) provided minor adaptations. 
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